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Abstract. This paper presents results concerning bifurcations of 2D 
piecewise-smooth vector fields. In particular, the generic unfoldings of 
codimension three fold-saddle singularities of Filippov systems, where a 
boundary-saddle and a fold coincide, are considered and the bifurcation 
diagrams exhibited. 

1. Introduction 

The general purpose of this paper is to study non-smooth vector fields 
(NSVF's for short), also called Filippov systems, Z represented by the fol- 
lowing three-parameter family of differential equations in M 2 : 

_, _ Xx= {a 1 (r)(x-X)la 2 (r)( X -Xf) if ^ °' 

where r = {inv,vis}, a\{inv) = — 1, a\{vis) = 1, a2(inv) = 1, ct2{vis) = 
and (A,/3 5/ u) G (-1,1) x (-^3/2, y%/2) x (-e ,l) with e > sufficiently 
small. 

In Figure [TJ (respectively, Figure [2]) we consider A = // = /? = and 
r = inv (respectively, r = vis) in Equation ([TJ. 

Call S x = {X 2 = 0} and S Y = \Y 2 = 0} where X x = X = (X\X 2 ) 
and = Y = (Y\Y 2 ). Denote T x ^ = (X X ,Y^) = (X,Y). In short 
our goal is to study the local dynamics of Z x p consisting of two smooth 
vector fields X\ and Y^ !( g in R 2 such that on one side of a smooth surface 
S = {y = 0} we take Z x g = X\ and on the other side Z^^p = Y^p- 
We mention two particular generic situations that occur in our system when 
(3 7^ 0. The first one is the fold- fold singularity (or two- fold singularity 
— see Definition [2]), studied in some detail in |Guardia et aL(2011)| and 
Kuznetsov et al.(2003)| ( among others), in which the trajectories of both 



X and Y have a quadratic tangency to S at a point p G E. The other 
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FIGURE 1. (Invisible) 
Fold-Saddle Singularity. 



FIGURE 2. (Visible) 
Fold-Saddle Singularity. 



situation is the occurrence of branches of canard cycles (see Definition [5]) 
which are typical minimal sets that appear in NSVF's. 

It is worth saying that the set fi r of all NSVF's Z = (X, Y) as below 
described (with some specified topology) is a differentiable manifold modeled 
in a Banach space (see Section [3] for details) . 

We now give a brief and rough overall description of the main features of 
this work: 

• In £l r all NSVF's presenting a generic fold-saddle singularity form a 
codimension-three submanifold W\ in such a way that any Z € W\ 
is locally equivalent (see Definition [8|) to the above described Zy„a 
with A = fj, = j3 = 0. 

• The bifurcation diagram of Z X fM /3 is exhibited for r = inv and r = 
vis. 

We emphasize that this paper is inserted in a larger 2D classification pro- 
gram where are included papers like [Teixeira(1977)||Teixeira(1991)[[Kuznetsov et aL(2003) 
Guardia et aL(2011)| . One of the goals of this program is to classify (via 
topological equivalence) typical singularities of NSVF's. For this purpose it 
is necessary to present generic unfoldings and give non-degeneracy condi- 
tions on the system in order to characterize the codimension of the singular- 
ity. The bifurcation diagram of the codimension-three singularity presented 
here includes that one exhibited in Guardia et aL(2011)| . The later claim is 
commented in the sequel. We finish this introduction presenting a physical 
model where a fold-saddle singularity can be found. 

Example 1. One of the reasons that second order linear equations with 
constant coefficients are worth studying is that they serve as mathemati- 
cal models of some important physical processes. Two important areas of 
application are in the fields of mechanical and electrical oscillations. For 
example, the motion of a mass on a vibrating spring, the torsional oscilla- 
tions of a shaft with a flywheel, the flow of electric current in a simple series 
circuit, and many other physical problems are all described by the solution 
of an equation of the form 



(2) 



ax + bx + cx = g(t, x). 
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Here we consider the external force g not depending on t but depending on 
x. For example consider 

(3) g(x) = Ax + 1 — sgn{x) with A > —. 

a 

Now if we call x = y then the equation ([2]) with g given by ([3]) became 

(4) x = y, y = (A - -) x - -y, if x > 
and 

(5) x = y, y = (A - -) x - -y + 2, if x < 0. 

V a) a 

(' 

System ([4]) has a saddle equilibrium at the origin because A > — . And 
system ([5]) has an invisible fold (see Definition {T^) at the origin. If g(x) = 
Ax — 1 + sgn(x) with A > — then it has a saddle equilibrium and a visible 

(X i — i 

fold. □ 

The paper is organized as follows: In Section [2] we present the main results 
of the paper. In Section[3]we give some basic concepts about NSVF's in order 
to setting the problem in Section HI The remaining sections are dedicated 
to prove the main results of the paper. 

2. Statement of the Main Results 
In what follows consider 

(6) notf) = 2 - (12/3/(-3 + 6/3 + ^9 - 12/3 2 )). 



Theorem 1. Assume t = inv and /U = /Uo(/3) in Equation {ip. Then its 
bifurcation diagram in the (A, ft)— plane contains 19 distinct phase portraits 
(see Fiaure [W\) . 



First of all observe that if /3 > and A = -1/2 + ^9 - 12/3 2 /6 in Theo- 
rem 1, then the NSVF has a homoclinic loop surrounding a non hyperbolic 
singularity. So, it is easy to see that the cases covered by Theorem 1 do not 
represent the full unfolding of the (invisible) fold-saddle singularity and the 
next two theorems become necessary. 

Theorem 2. Assume t = inv and Ho(/3) < fx < 1 in Equation fip. Then its 
bifurcation diagram in the (A, ft)— plane contains 21 distinct phase portraits 
(see Fiaure [2l\) . 

Theorem 3. Assume r = inv and —Eq < fi < /j-o(/3) in Equation |7p. 
Then its bifurcation diagram in the (A, /?)— plane contains 21 distinct phase 
portraits (see Figure [21]) . 
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Remark 1. The bifurcation diagrams exhibited in Theorems 2 and 3 present 
a homoclinic loop surrounding a hyperbolic singularity. Observe that, under 
the conditions of Theorem 2 (respectively, Theorem 3) this singularity is an 
attractor (respectively, repeller) as illustrated in Figure^ cases (a) and (b). 
Moreover, when the parameter ^ varies from —Eq to 1 there is an element 
Ho(f3) £ (— £0) 1)> given by Equation |6p, such that Z^^/pyp presents a like 
Hopf bifurcation phenomenon as illustrated in Figure This phenomenon 
is fully treated in |Guardia et al. (2011)] and Kuznetsov et al. (2003) | and it 
is not covered by Theorems 1, 2 and 3. 






(a) 



(b) 



(c) 



FIGURE 3. A like Hopf bifurcation. The singularities in (a) and 
(c) (that appear in Theorems 2 and 3 respectively) are called 
S— attractor and S— repeller throughout this paper (see Figure [5] 
cases (b) and (c)). In (b) we observe a non hyperbolic singularity 
present in Theorem 1. 



Remark 2. In Theorems 2 and 3 beyond the cases (a) and (c) in Figure 
it also appears another relevant phenomenon that is the presence of a 
fold-fold singularity. Note that in Theorem 1 it occurs simultaneously the 
occurrence of a loop (a global phenomenon) and a fold-fold singularity (a 
local phenomenon). 

Theorem 4. Assume r = vis in Equation or equivalently, take r = vis 
and [i = in Equation |7]). Then its bifurcation diagram in the (A, /3)— plane 
contains 13 distinct phase portraits (see Figure {27\). 

The cases covered by Theorem 4 do not represent the full unfolding of the 
(visible) fold-saddle singularity. So, the next two theorems are necessary. 
Each one of them describes a distinct generic codimension two singularity. 

Theorem 5. Assume r = vis and < fj, < 1 in Equation (Op. Then its 
bifurcation diagram in the (A, /?)— plane contains 13 distinct phase portraits 
(see Figure{27\). 
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Theorem 6. Assume t = vis and £q < fi < in Equation Then its 

bifurcation diagram in the (A, /3)— plane contains 13 distinct phase portraits 
(see Figure[27\). 

Remark 3. In Theorems 5 and 6 one observes the birth of a singularity on 
£. This singularity behaves like a saddle (see Figure\£&-(a)) and is known as 
£— saddle. This phenomenon does not occur under the conditions presented 
in Theorem 4- 

3. Basic Theory about NSVF 

Let K C R 2 be a compact set such that dK is a smooth curve. Con- 
sider SCii" given by £ = / _1 (0), where / : K — > R is a smooth function 
having G R as a regular value (i.e. V/(p) 7^ 0, for any p G 7 _1 (0)) 
such that dK n £ = or dK rti £. Clearly the switching manifold £ 
is the separating boundary of the regions £ + = {q £ K \ f(q) > 0} and 
XL = {q £ K I /(q) < 0}. We can assume that E is represented, locally 
around a point q = (x, y), by the function f(x, y) = y. 

Designate by x r tli e space of C r — vector fields on K endowed with the 
C r — topology with r > 1 large enough for our purposes. Call O r = Q r (i^, /) 
the space of vector fields Z : K — > R 2 such that 

for 0,y)GE+, 
YO,y), for (a?,y)e£_, 

where X = (X 1 ,^ 2 ) and Y = (Y\Y 2 ) are in X T '■ We write Z = (X,Y), 
which we will accept to be multivalued in points of E. The trajectories of 
Z are solutions of q = Z(q), which has, in general, discontinuous right- 
hand side. The basic results of differential equations, in this context, were 
stated by Filip pov in [Filippov(1988)] . Related theories can be found in 
[Kozlova(1984)| |Teixeira(1991)] among others. 



Z(x,y) 



3.1. Orbits, trajectories and singularities of NSVF's. In what follows 
we will use the notation 

X.f(p) = (Vf(p),X(p)) and X\f(p) = (V^" 1 ./)^), X(p)) , i > 2 

where (., .) is the usual inner product in IR 2 . 

Following the Filippov rule, we distinguish the following regions on the 
discontinuity set E : 

• Crossing region: E c = {p G E | (X.f(p))(Y.f(p)) > 0}. 

• Escaping region: E e = {p G E | X.f(p) > and Y.f(p) < 0}. 

• Sliding region: E s = {p £ E | X.f(p) < and Y.f(p) > 0}. 

Consider Z = (X, Y) G VL r and p G E e U E s . The sliding vector field Z s 
associated to Z at p is the convex combination of X(p) and Y(p) tangent to 
E at p (see Figure H} . 
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FIGURE 4. Filippov's convention. 



We say that q G X is a X— regular point if 
. (X.f(q))(y.f(q)) > or 

• (X.f(q))(Y.f(q)) < and Z s (g) / (that is q G X e U X s and it is 
not an equilibrium point of 

The points of X which are not X— regular are called T,— singular. We dis- 
tinguish two subsets in the set of X— singular points: X* and X p . Any q G X p 
is called a pseudo equilibrium of Z and it is characterized by -Z s (o) = 0. Any 
g G X* is called a tangential singularity and is characterized by Z T '(q) ^ 
and (X/(g))(Y./(g))=0. 

A pseudo equilibrium g £ 57 is a X— saddle provided that one of the 
following condition is satisfied: (i) 5 6 E e and q is an attr actor for or 
(ii) 5 £ S s and q is a repeller for Z s . A pseudo equilibrium q G X p is a 
T,— repeller (resp. X— attractor) provided g 6 S e (resp. g G X s ) and g is a 
repeller (resp. attractor) equilibrium point for Z s (see Figure [5]). 




(a) (b) (c) 

FIGURE 5. Pseudo equilibria: (a) £— saddle, (b) £— attractor and 
(c) repeller. 



Definition 1. We say that po G X* is a "E—fold point of X G x r if 

X.f(po) = but X 2 .f(po) 7^ 0. Moreover, po G X is a visible (respectively 
invisible) X— /o/d point of X if X.f(po) = and X 2 ./(po) > (respectively 
X 2 .f( Po ) < 0). 
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Definition 2. Let Z = (X, Y) G Q r . We say that p G E* is a fold- fold 
singularity of Z if p is a H—fold point for both X and Y . 

Definition 3. Let Z = (X,Y) G Q r . We say that q G E is a fold-saddle 
singularity of Z if q is a T,— fold point of X and a saddle equilibrium ofY 
(in this case q is called a boundary-saddle ofY). 

The following construction is presented in |Teixeira(2008)| . Let Z = 
(X, Y) G f2 r such that T is an invisible E— fold point of X. From Implicit 
Function Theorem, for each p G E in a neighborhood Vt of T we derive 
that there exists a unique t(p) such that the orbit— solution t \— > tpx(t,p) 
of X through p meets E at a point p = ipx(t(p),p)- Define the map 
Ix ■ Vt H S — > Vt H E by jx(p) = P- This map is a C— diffeomorphism 
and satisfies: 7^ = Ld. Analogously, when T is an invisible E— fold point of 
Y we define the map 7y : Vf n E — > Vf n E associated to Y which satisfies 
7y = Ld. We define now the first return map associated to Z = (X, Y): 

Definition 4. The first return map ipz '■ T — > T is defined by the com- 
position tpz = 7y o 7x when both and 7y are well defined in T C E. 

Remark 4. ipz is an area-preserving map. 

Definition 5. A curve T is a canard cycle of Z = (X, Y) G O r if it is 

closed and composed by orbit-arcs of at least two of the vector fields X\-£ + , 
Y\-£_ and Z s . We say that T is hyperbolic if ip' z (p) ^ 1, where ipz is the 
first return map defined on a segment T with p G T iti 7. 

Definition 6. Consider Z G Sl r . ^4 closed path A is a E— graph if it is a 

union of equilibria, pseudo equilibria, tangential singularities of Z and orbit- 
arcs of Z joining these points in such a way that A n E 7^ 0. 

Definition 7. Consider Z £ Q r . A point q G E is a E— center if there is 
a neighborhood U of q filled up with a one-parameter family of canard 
cycles of Z in such a way that ju n E C S c (see Figure^). 





Figure 

6 . Canard 
cycle. 



Figure 7. E-graph 

surrounding a E— center. 
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3.2. Structural Stability on £l r . Bifurcation theory describes how con- 
tinuous variations of parameter values in a dynamical system can, through 
topological changes, cause the phase portrait to change suddenly. In this 
paper we focus on certain structurally unstable NSVF's within a generic 
context. In Andronov Sz Pontryagin(1937)| the concept of k^— order struc- 



tural stability is presented; in a local approach such setting gives rise to the 
notion of a codimension k singularity. Now we present the concept of equiv- 
alence which will guide us for all the paper. 

Definition 8. Two NSVF's Z, Z G n r (K, f) defined in open sets U, U C K 
with switching manifolds T, C U and £ C U respectively are XI— equivalent 
if there exists an orientation preserving homeomorphism h : U —> U which 
sends £ to S and sends orbits of Z (respectively Z") to orbits of Z (respec- 
tively Z^). From this definition the concept of local structural stability 
in Q r is naturally obtained. 

As we said in Section Q] our paper is a generalization of some papers that 
unfold typical singularities of NSVF's. Below we present the program used 
in the literature and in our paper to exhibit the diagram bifurcation of a 
singularity of NSVF's. 

Let Z G J7 r andp G S. Following the approach in |Sotomayor(1974j| (and 
also exposed in Guardia et q/.(2011)| , Subsection 3.1, pg 1982), we get that: 



By Theorem 3.5 of [Guardia et aZ.(2011)] we already know the char- 
acterization of the set $q = {Z G 0, r | Z is locally structurally stable at p} . 
In fact, $q is open and dense in £l r . So, is the codimension zero 
local bifurcation set. 

Let S7i = r2 r \$Q and = {Z G fii | Z is locally structurally stable at p relative to fii}. 
The set is the codimension 1 local bifurcation set in f2 r . The 
characterization was given in [Kuznetsov et q/.(2003)]. 

In addition, if Zq G $f , it is also known (see |Teixeira(1979)| ) that 
there exists a neighborhood U of Zq in Q r such that: 

o There exists a C r — function L : U — > R, satisfying L(Zq) = 
and DLz , the differential of L at Zq, is surjective. Moreover, 

o Consider now all the embeddings : (— e, e) x U — > U 
transversal to at some Z G $J and such that 0(0, Zq) = Z. We 
refer to such as an unfolding of Zq. 

We consider now the set 0,2 = an d similar objects <3?2 (the 

set of codimension 2 singularities) and families of objects L : U — > 
M?, with surjective derivative at Zq and embeddings : (— e, e) x 

(-CC) xu^u. 

In this way we can get sequences fi^ and $| in fi r , that establish a 
program to characterize all codimension k singularities. 
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Definition 9. Let V(0,R fc ) and 5(0, R l ) be neighborhoods of in R k and 
W respectively and let U be a neighborhood of Zq in il r . We say that two 
unfoldings 6 : V(0, R fc ) x U — > U and E : 5(0, M. 1 ) xU — >• U are equivalent 
if there is a homomorphism A : V(0, R fc ) — > 5(0, M. 1 ) such that A(X) = fi and 
for each Z G the vector fields Q(X,Z) and z,(A(X),Z) are S— equivalent 
according to Definition^ Moreover, we say that an unfolding @(Ao, •) is a 
generic unfolding if there is a neighborhood W(G(Ao,-)) °/0(^o>-) such 
that any unfolding 0(A, .) G W(G(Ao, •)) i> s equivalent to 0(Ao, •)• 

Remark 5. In Definition^ it is important to say that the homomorphism 
A does not vary, necessarily, continuously with respect to the parameter 
A G V(0,M fc ). 

3.3. The Direction Function H. Here we introduce a function that will 
be very usefull in the sequel. 

In (A, B) C E e US s , consider the point C = (Ci, C2), the vectors X(C) = 
{D 1 ,D 2 ) and Y(C) = (E 1 ,E 2 ) (as illustrated in Figure ED- The straight 
segment passing through C + X(C) and C + Y(C) meets S in a point p(C). 
We define the C r — map 

p: (A, B) — > E 
z 1 — > p(z). 

Since £ is the x— axis, we have that C = (Ci, 0) and p(C) 6ix {0} can be 
identified with points in R. According with this identification, the direction 
function on E is defined by 

H : (A,B) — > R 

z 1 — >• p(z) — z. 



X 



A 
Y 




C + Y(C) 



B 

C + X(C) 



FIGURE 8. Direction function. 



We obtain that H is a C r — map and 

• if H(C) < then the orientation of Z s in a small neighborhood of 
C is from B to A; 

• if H{C) = then C G £P; 

• if H(C) > then the orientation of Z s in a small neighborhood of 
C is from A to £?. 



10 



C.A. BUZZI, T. DE CARVALHO AND M.A. TEIXEIRA 



Simple calculations show that p(C\) = E2 ^ Dl+C ^ 2 _ D ]^ El+Cl ^ and conse- 
quently, 

(7\ u(r \ - E2Dl Z D2El 

(?) H{Cl) ~ E 2 -D 2 ■ 

Remark 6. If X.f(p) = and Y.f(p) ^ then, in a neighborhood V p of p in 
£, the direction function H has the same sign of D\, where X(p) = (D\, D 2 ). 
In fact, X.f(p) = and Y.f(p) 7^ are equivalent to D 2 = and E 2 ^ in 

j*7|). So, lim H(pi) = D\, where ko 7^ and p = (pi,p 2 ). 

(D 2 ,E 2 )^(0,k ) 



4. Setting the problem 

Let T+ = {X\ s+ with X £ xl (respectively, T_ = {X| s _ with X £ 
X r })- This means that r + (respectively, T_) is identified with % r . 

Let C r + be the set of all elements X € r + having a S— fold point, 
is an open set in r + (see [Teixeira(1977)1 ). We may consider f(x,y)=y 
and the following generic normal forms Xo(x,y) = {a\,(3\x) with a± = ±1 
and Pi = ±1 (see [Vishik(1972)1 , Theorem 2). 

Let r|_ C T_ be the set of all elements Y G T_ presenting a hyperbolic 
saddle equilibrium Sy on E (called boundary saddle in the literature — as 
a reference, in |Roy fc Roy (2008)] , are exhibited some border collisions in 
three-dimension NSVF's) and such that the eigenspaces of DYq(Sy ) are 
transverse to £ at Sy ■ is a codimension one submanifold of T_ . Note 
that, has a C r — structure (for more details about this construction see 
Dumortier (1978)] ). In fact, since Yq has a hyperbolic boundary saddle 
Sy and the eigenspaces of DYq(Sy ) are transverse to £ at SV , then there 
exists a neighborhood V(io) of Yq in x r such that all Y G V(lb) have 
a hyperbolic boundary saddle Sy = (x*,y*) with the same properties of 
the eigenspaces. Moreover, the correspondence I : V(Yq) — > M 2 , given 
by L(Y) = Sy, is C r at Y . Define II : V(Y ) M, given by U(Y) = 
(f o L){Y). We say that two vector fields Y,Y € defined in open 

sets U and U, respectively, are C°— orbitally equivalent if there exists an 
orientation preserving homeomorphism h : U — >■ U that sends orbits of 
Y to orbits of Y. From |Teixeira(1977)| we know that any 7 £ r| is 

generically C°— orbitally equivalent to its linear part by a S— preserving 
homeomorphism. And the linear saddle with eigenspaces transverse to the 
x— axis has the generic normal forms Yo(x,y) = (a 2 y,a 2 x) with a 2 = ±1. 
So it is easy to see that the generic unfolding of the singularity is given by 

= (02(11 + P), cn 2 x) where /3 G R. 

Moreover: 

• There exists a C r — function II : V(Yq) — > R, such that -DIIy is 
surjective. 

• The correspondence Y 1— > Sy is C r , where Sy is a saddle point of Y. 
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• If II(Y) < then S Y G E_ 

• If U(Y) = then S Y G E. 

• If II(Y) > then G E+ 



In this paper we are concerned with the bifurcation diagram of systems 
Zq = (Xq,Yq) in Q r such that Xq G ^0 G T|_ and po = Sy G E. 

The fold-saddle singularity po = Sy is illustrated in Figures Q] and [2] — the 
dotted lines in these and later figures represent the points where X.f = 
and Y.f = 0. 

Let p = (0,0) be a fold-saddle singularity of Z = (X, Y). We denote the 
set of all NSVF Z = (X, Y) such that X £ r£ + and Y G rg by r F ~ 5 . We 

endow % r x \ r (consequently £l r and T^ 5 ) with the product topology. Let 
Z = (Xq, Yq) E r F " 5 . Observe that is the unique singularity of Xq around 
a neighborhood Wq of the origin in M 2 . So, there exists a neighborhood C/o 
of Zq in such that for any Z = (X, Y) £ C/o we may find a E— fold 
point pz = (kz, 0) G Wq such that it is the unique singularity of X in Wo. 
Moreover the correspondence Z ^ pz is C r . 

In the same way, for any Z = (X, Y) G Uq we find a C r — correspondence 
B : C/o — J- M? where B(Z) = sz = (az,bz) is the (unique) equilibrium 
(saddle) of Y in C/o- We are assuming that the eigenspaces of DY Sz (qz) are 
transverse to E at sz- We have to distinguish the cases: (i) bz < 0, (ii) 
bz = and (m) bz > 0. Observe that when bz < (resp. &z > 0) there 
is associated to Z an invisible (resp. visible) E— fold point of Y given by 
qy = (cz,0) G Wo- Moreover limi >z ^o c z = a>z- 

Define F(Z) = (kz — az,bz)- Knowing that sz is a hyperbolic equilib- 
rium point of Z and is a regular value of F, it is not hard to prove (see 
[Teixeira(1979)l ) that: 

• The derivative DF : Uq — > M 2 is surjective and 

• -F -1 (0) = 0,2 is a codimension two submanifold of f2 r . 

Therefore this fold-saddle singularity occurs generically in two-parameter 
families of vector fields in Q r . 

4.1. Normal Form. We start this section with the following model: 



where r = {inv,vis}, a\{inv) = —1, a\{vis) = 1, p\ = ±1 and k\ = ±1. 

The next lemma provides explicitly the equivalence between Z S and 
the model ([8]). We present an outline of proof of the previous lemma in 
Section [5j 




Lemma 1. If Z € Q2 then Z is E— equivalent to Z T given by (Gjj. 
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Note that there exists a NSVF Zeff nearby Z inv , given by ©, such 
that Z presents a S— center (see Figure [7j). In fact, this suggests that the 
unfolding of (JSj) has infinite co dimension. At this point it seems natural to 
propose the following conjecture. 

Conjecture: For any neighborhood W C O r of Z mv (given by $Bj)), and 

for any integer k > there exists Z S W such that the codimension of Z is k. 

So, based on this conjecture, we have to sharpen our normal normal. In 
fact, in order to get low codimension bifurcation we have to impose some 
generic assumptions. 

Without loss of generality, throughout the rest of this paper we consider 
pi = 1 and k\ = — 1 at the model ([8]). The other choices of the parameters 
are treated analogously. When r = inv we add the extra generic assumption 
X$.f(p) ^ on the S-fold po int p of Z = (X ,Y ) G T F ~ S . By means 
of Theorem 2 in [Vishik(1972)] , we may conclude that around the invisible 
S— fold point the vector field Xq can be expressed as Xq = (1, — x + a\x 2 ) 
and f(x, y) = y, where a\ / 0. We say that Xq is contractive (respectively, 
expansive) at p if a\ < (respectively a\ > 0). 

According to the previous discussion, we will consider Z^^Zq 113 E QJ' 
written in the following forms: 



Note that X 1 ™ presents an invisible expansive S— fold point in its phase 
portrait and Xq is presents a visible one. 

4.2. Unfoldings. The main question of this paper is to exhibit the bifur- 
cation diagram of Z$ with either r = inv or r = vis. For this reason we 
consider unfoldings of the normal forms ([9]) and ()10p . We obtain that: 



I- There is an imbedding F r : M 2 ,0 -> x r ,^o such that F o( x ^) = z \. 
is expressed by: 



(9) 




(10) 




(11) 
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where r = {inv,vis}, (A,/3) G (—1,1) x (— \/3/2, y/3/2), ai(inv) = — 1, 
ai(ws) = 1, a2{inv) = 1 and ct2(vis) = 0. Moreover, the two-parameter 
family given by (jlip is transversal to We stress that in Guardia et a/.(2011)| 
the unfolding of the case r = inv is done. Nevertheless we observe that there 
are some typical topological types nearby Zq that do not appear in the bi- 
furcation diagram of Zj^ g. For example, when r = inv the configurations in 
Figure [9] cases (a) and (b) are excluded and when r = vis the configuration 
in Figure 0- (c) also is excluded. 




(a) (b) (c) 

FIGURE 9. Cases do not covered in |Guardia et aZ.(2011)| . In 
(a) and (b) (respectively, in (c)) the fold singularity is invisible 
(respectively, visible) . 



II- In order to consider a more general situation which includes the above 
mentioned cases we add an auxiliary parameter u. As a result the model in 
Equation [T] is obtained. 

We stress that the configuration illustrated in Figure [H— (a) plays a very 
important role in our analysis. In this resonant configuration we note, si- 
multaneously, a fold-fold singularity (which is a local phenomenon) and a 
E— graph (loop) passing through the saddle equilibrium (which is a global 
phenomenon). Only the bifurcation of these two unstable configurations al- 
ready represents a relevant development (the fold-fold singularity was stud- 
ied recently in |Guardia et a?.(2011)| and the non-smooth loop bifurcation, 
as far as we know, was not studied until the present work). In fact, this 
configuration is reached in ([T]), taking n = n (/3). 
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In what follows, in order to simplify the calculations, we take fx = a + 1 
in ([I]) and obtain the following expression: 



fl2) 



Y 



a,/3 



ai(r)(x - A) + a 2 (r)(x - A) 
k ^W^(y + /3) 



if y > o, 
if y < o, 



where r = {mu, uis}, ai(inu) = —1, ai(vis) = 1, a2(inv) = 1, a^vis) = 
and (A,/3,a) G (-1, 1) x (- y/3/2, y/3/2) x (-l-e , 1) with e > sufficiently 
small. Since [J>o(/3) is given by ©, we obtain that 



(13) 



a 



{P) = 1 - (12/3/(-3 + 6/3 + V9-12^ 2 )). 



When it does not produce confusion, in order to simplify the notation we 
use Z = (X,Y) or Z x>a ,p = (X,Y) instead Z\ a ^ = (X x ,Y a>/3 ). 



4.2.1. Geometrical Analysis of the Normal Form \12}) . Given Z = (X,Y), 
we describe some properties of both X = X\ and Y = Y a> 0. 

The real number A measures how the S— fold point al = (^1,^2) = (A,0) 
of X is translated away from the origin. More specifically, if A < then d 
is translated to the left hand side and if A > then d is translated to the 
right hand side. 

Some calculations show that the curve Y.f = is given by y = — /3. 

So the points of this curve are equidistant from the separatrices when a = 
— 1. It becomes closer to the stable separatrix of the saddle equilibrium 
S = = (s±,S2) when a E (—1,0). It becomes closer to the unstable 
separatrix of S when a G (—1+Eq,—1). Moreover, the smooth vector field Y 
has distinct types of contact with £ according to the particular deformation 
considered. In this way, we have to consider the following behaviors: 

• Y~ : In this case j3 < 0. So S is translated to the y— direction with 
y > (and S is not visible for Z). It has a visible S— fold point 

e = e a> p = ( ei ,e 2 ) = ({gf/3,o) = (ei,0) (see Figure HDD- 

• Y° : In this case (3 = 0. So S is not translated (see Figure [1]). 

• Y + : In this case (3 > 0. So S is translated to the y— direction 
with y < 0. It has an invisible S— fold point i = i a ^ = (11,12) = 

( (i-°) ft' 0) ' Moreover, we distinguish two points: h = hp = (hi, /12) 
= (— j3, 0) which is the intersection between the unstable separatrix 
with E and j = jp = (j 1,^*2) = (AO) which is the intersection 
between the stable separatrix with S (see Figure [TT1) . 

In Figure [TT] we distinguish the arcs of trajectory a\ joining S to h and 
a 2 joining j to S. 
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Figure 10. Case Figure 11. Case 

Y~: perturbation of Y + : perturbation of 

a boundary saddle a boundary saddle 

generating a visible generating an invisible 

fold. fold. 



5. Proof of Lemma 1 
Now, we show how we can construct the homeomorphism in Lemma [TJ 



Outline of Proof of Lemma 1. Here we construct a £— preserving homeo- 
morphism h that sends orbits of Z = (X,Y) € f^j defined in a suffi- 
ciently small neighborhood Uq of the fold-cusp singularity of Z, to orbits of 
Z = (X , Y), defined in a sufficiently small neighborhood Uq of the fold-cusp 
singularity of Z, where Z = Z inv is given by (jSJ) with pi = 1 and k\ = — 1 
(the other cases are treated analogously). Consider Aq an arbitrary point 
of the stable separatrix of the saddle point S of Y (see Figure [T2|) . Let Ti 
be a transversal section of Y at Aq. The section T\ also is transversal to Y 
and it crosses the stable separatrix of the saddle point S of Y at -Bo- Let 
A\ € T\ be a point on the left of Aq. The trajectory of Y passing through 
A\ crosses £ at A2. In the same way, the trajectory of Y passing through 
B\ crosses £ at B<i- The trajectory of X that passes through A2 crosses 
£ in a point A3. The trajectory of X that passes through B2 crosses £ in 
-B3. Consider A^ an arbitrary point of the unstable separatrix of S. Let 
T2 be a transversal section of Y passing through A±. The section T2 also 
is transversal to Y and it crosses the unstable separatrix of S at the point 
B4. The trajectory of Y passing through A3 crosses T2 in a point A§. In 
the same way, the trajectory of Y passing through B3 crosses T2 at B$. Let 
Aq G Ti be a point at the right of Aq. The trajectory of Y passing through 
Aq crosses T2 at A7. The trajectory of Y passing through Aq crosses T2 at 
B-j. The homeomorphism h sends T\ to T±, the arc of trajectory 71 = A\ A5 
to the arc of trajectory 71 = A\ B§ and the arc of trajectory 72 = Aq A7 to 
the arc of trajectory 72 = AqBj. Now we can extend continuously h to the 
interior of the region bounded by T\ U 71 U T2 U 72. In this way, there exists 
a E— preserving homeomorphism h that sends orbits of Z to orbits of Z. 
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When Z = Z ms is given by (JSj) , the first coordinate of X is equal to 1 and 
k\ = —1. We build the same above construction until the appearance of A2 
and I?2- Now, consider C4 an arbitrary point of the stable separatrix of the 
X— fold point F of X (see Figure fl3|) . Let T3 be a transversal section to 
X at C4. The section T3 is also transversal to X and it crosses the stable 
separatrix of the S— fold point F of X at the point D4. The trajectory 
of X passing through Ai crosses T3 at a point C3. In the same way, the 
trajectory of X passing through B2 crosses T3 at a point D3. Consider C5 
an arbitrary point of the unstable separatrix of F. Let T4 be a transversal 
section to X at C5. The section T4 also is transversal to X and it crosses 
the unstable separatrix of F at D5. Let Cq E T3 be a point at the right 
of C4. The trajectory of X passing through Cq crosses T4 at C7. In the 
same way, let Z?6 G T3 be a point at the right of D4. The trajectory of X 
passing through Dq crosses T4 at Dj. Let Cg S T4 be a point at the right of 
C5. The trajectory of X passing through Cg crosses £ at ^3. In the same 
way, let -Dg S T4 be a point at the right of -D5. The trajectory of X passing 
through D$ crosses S at S3. Now, it is enough to repeat the construction 
made in the previous case. The homeomorphism h sends T\ to 7\, the arc 
of trajectory 71 = A\ C3 to the arc of trajectory 71 = A\ D3, T3 to T3, the 
arc of trajectory 72 = Cq C-j to the arc of trajectory 72 = Dj, T4 to T4, 
the arc of trajectory 73 = Cs A§ to the arc of trajectory 73 = D$ B$, T2 to 
T2 and the arc of trajectory 74 = Aq A-j to the arc of trajectory 74 = Aq Bj. 
Now we can extend continuously h to the interior of the region bounded by 
T\ U 71 U T3 U 72 U T4 U 73 U T2 U 74. In this way, there exists a S— preserving 
homeomorphism h that sends orbits of Z to orbits of Z . □ 



ON THREE-PARAMETER FAMILIES OF FILIPPOV SYSTEMS - THE FOLD-SADDLE SINGULARITY 




FIGURE 13. Construction of a C°— equivalence: Case r = vis. 



6. Proof of Theorem 1 

Proof of Theorem 1. In Cases li, 2\ and 3i we assume that Y presents the 
behavior Y~ . In Cases 4i, 5i and 6i we assume that Y presents the behavior 
Y°. In these cases canard cycles do not arise. 

o Cases (li) d\ < e±, (2i) d\ = e\ and (3i) d\ > e±: The points of £ 
outside the interval (d±,e±) (or (e±,di)) belong to £ c . The points inside this 
interval, when it is not degenerate, belong to £ s in Case li and to £2 i n 
Case 3i. In both cases H(z) > for all z£S e U S s . See Figure [T4l 



A<0 A=0 A>0 




o Cases (4i) d\ < s\, (5i) di = si and (61) d\ > s±: The points of £ 
outside the interval (di, si) (or (si, di)) belong to S c . The points inside this 
interval, when it is not degenerate, belong to T, s in Case 4\ and to £2 hi 
Case 61. In both cases H(z) > for all z E £ e U £ s . See Figure [HJ 

In Cases 7\ — 19i we assume that Y presents the behavior Y + . 

Observe that as the parameter A increases — assuming the values Lq, L\ 
and L2 below described — it appear orbit-arcs of X connecting the points h 
and i, h and j and i and j respectively. 

Remembering that a = 1 - (12/3/(-3 + 6/3 + ^9 - 12/3 2 )) the values of 
Lq, Li and L2 are: 
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A<0 A=0 A>0 



4i ; 




'4\ 





















FIGURE 15. Cases 4i, 5i and 6i. 

L = [-9 - 6/3 + ^9 - 12/3 2 + 15 + ^9 - 12/5 2 - 2/3(-3 + 2/3 + v^^l^)]/^ 
Li = -1/2 + ^9-12/32/6 

L 2 = [-9 + 6/3 + ^9- 12/3 2 + 15 + ^9 - 12/3 2 + 2/3(-3 + 2/3 + ^9^12^)]/ 12. 
o Cases (7i) A < -/3, (8i) A = -0, (9i) -/3 < A < L , (10i) A = L 
and (Hi) Lq < A < L\: The points of E outside the interval (di, i\) belong 
to E c . The points inside this interval belong to E s . The direction function 
if assumes positive values in a neighborhood of di, negative values in a 
neighborhood of i\ and there exists only one value P = P\ a r such that 
ff(P) = 0. So, by ([7]), the E-attractor P = (P,0), nearby (0,0), is the 
unique pseudo equilibrium of Z. In these cases canard cycles do not arise. 
See Figure [HJ 




o Case (12i) A = L%: Since A = L± there is an orbit-arc 7^ of X con- 
necting the points h and j. It generates a E— graph V = 75^ U E e U 5 U E c 
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of kind I. Moreover, since a = oq, where ao is given by f)13f) . there exists 
a non generic tangential singularity at the point d = i. So, the points of 
Y,/{d} belong to S c . As the S— fold point of X is expansive, a direct calculus 
shows that the First Return Map r\ : (h, d) — > (h, d) has derivative bigger 
than 1. As consequence, T is a repeller for the trajectories inside it, d = i 
behaves itself like an attractor (weak focus) and canard cycles do not arise. 
See Figure [I7j 




o Case (13i) L\ < A < L3; The meaning of L3 will be given below in 
this case. The points of £ outside the interval (ii, d\) belong to S c and the 
points inside this interval belong to S e . The direction function H assumes 
positive values in a neighborhood of d\, negative values in a neighborhood 
of i\ and there exists a unique value P = P\,a,p such that H(P) = 0. So 
P = (P,0) is a £— repeller. When A is a bit bigger than L\, the First 
Return Map r\ has two fixed points, i.e., Z has two canard cycles. One 
of them, called T\, is born from the bifurcation of the S— graph T of the 
previous case and the other one, called T2, is born from the bifurcation of the 
non generic tangential singularity presented in the previous case. Both of 
them are canard cycles of kind I. Moreover, we obtain that T± is a hyperbolic 
repeller canard cycle and T2 is a hyperbolic attractor canard cycle. Note 
that, as A increases, T± becomes smaller and T2 becomes bigger. When A 
assumes the limit value L3, one of them collides to the other. See Figure [T71 

o Case (14i) A = L3; The distribution of the connected components of T, 
and the behavior of H are the same as Case 13i. Since A = L3, as described 
in the previous case, there exists a non hyperbolic canard cycle T of kind I 
which is an attractor for the trajectories inside it and is a repeller for the 
trajectories outside it. See Figure [UJ 

o Cases (15i) L 3 < A < L 2 , (I61) A = L 2 , (17i) L 2 <\</3, (I81) A = /3 
and (19i) A > (3: The points of S outside the interval belong to 

S c and the points inside this interval belong to E e . The direction function 



20 



C.A. BUZZI, T. DE CARVALHO AND M.A. TEIXEIRA 



H assumes positive values in a neighborhood of d±, negative values in a 
neighborhood of i\ and there exists a unique value P such that H(P) = 0. 
So, by 0, the E— repeller P = (P, 0), nearby (0,0), is the unique pseudo 
equilibrium of Z. In these cases canard cycles do not arise. See Figure [T3 

L 3 <\<L 2 X = L 2 L 2 <\<(3 \ = /3 X > 




Figure 18. Cases 15i - 19i. 




I 

FIGURE 19. Bifurcation Diagram of Theorem 1. 



The bifurcation diagram is illustrated in Figure [T9j □ 

Remark 7. In Cases Hi and 15i the ST -bifurcations (as described in 
Guardia et aZ.(2011)| ; arise. In fact, note that the trajectory passing through 
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h can make more and more turns around P. This fact characterizes a global 
bifurcation also reached in other cases. 

7. Proof of Theorem 2 

Proof of Theorem 2. In Cases I2, 22 and 32 we assume that Y presents the 
behavior Y~ . In Cases 42, 52 and 62 we assume that Y presents the behavior 
Y°. In Cases 7 2 - 21 2 we assume that Y presents the behavior Y + . 

o Cases (1 2 ) di < e\, (2 2 ) d\ = e\, (3 2 ) d\ > e\, (4 2 ) d\ < s\, (5 2 ) 
d\ = s\ and (62) d\ > s±: The analysis of these cases are done in a similar 
way as the cases li, 2i, 3i, 4i, 5i and 61. 

Observe that as the parameter A increases — assuming the values Mq, Mi 
and M 2 below described — it appear orbit-arcs of X connecting the points 
h and i, h and j and i and j respectively. The values of Mq, M\ and M2 
are: 

M = (-3 - 3a(-2 + a + 2(-l + a)/3) + 

+ V9(-l + a) 4 - 1 2(-1 + a) 2 ^)/(6(-l + a) 2 ), 
Mi = -l/2 + ^9-12/3 2 /6 
M 2 = (-3 + 6/3 - 3a(-2 + a + 2/3)+ 

+ ^/9(-l + a) 4 - 12(-1 + a) 2 a 2 /3 2 )/(6(-l + a) 2 ), 
o Cases (7 2 ) A < -(3, (8 2 ) A = -0, (9 2 ) -/3 < A < Af , (10 2 ) A = M 
and (H2) Mq < A < Mi: Analogous to Cases 7i — Hi changing Lq by Mq 
and L\ by Mi. 

o Case (122) A = Mi: The points of £ outside the interval {d\,%\) belong 
to S c and the points inside this interval belong to S s . The direction function 
H assumes positive values in a neighborhood of di, negative values in a 
neighborhood of i\ (see Remark[6]) and there exists a unique value P = P\, a ,/3 
such that H{P) = 0. So P = (P, 0) is a S-attractor. Since A = Mi, there is 
an orbit-arc of X connecting the points h and j. It generates a S— graph 
T = 7^ U S e U 5U S c of kind I. Since a > ao, where ao is given by ([TBI) , it is 
straightforward to show that the First Return Map defined in the interval 
(hi,di) C S do not have fixed points. By consequence, T is a repeller for 
the trajectories inside it and canard cycles do not arise. See Figure 1201 

o Case (132) Mi < A < i\: The distribution of the connected components 
of X and the behavior of H are the same as Case 122- Since Mi < A < i\, 
there is an orbit-arc 7^ of X connecting j to a point k = (fci,0) £ S, 
where k\ € (h±,di), for negative time. Also there is an orbit-arc 7^ of Y 
connecting A; to a point I = (h,0) € E, where Zi £ for negative time. 

Repeating this argument, we can find an increasing sequence We 
can prove that there is an interval / C (k, d) such that 7/ = (<£>y o ipx)' < 1 
on I. As P is a X— attractor, there is an interval J C d) such that 
rj' > 1 on J. Moreover, we can prove that 77 has a unique fixed point 
Q G (k,d). As consequence, there passes through Q a repeller canard cycle 
r of kind I. See Figure [20l This canard cycle is born from the bifurcation of 
the 2— graph present in Case 122- The expression of 77 is too large, so the 
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A = Mi Mi < A < h A 




Figure 20. Cases 12 2 , 13 2 and 14 2 



general case will be omitted. For the particular case when a 
and A = -1/2 + 11^6/60, the map r) is given by 



■1,0 = 1/2 



T](x) 



3,3 
4 + 2 



2 + lOy^ + 2^ 



1 



2 r lOv^ 



2x 3 + 2 



1 — L_ -U_ 

2 ^ lOv^ 



+ 2x 



A straightforward calculation shows that the unique fixed point of this par- 
ticular rj occurs when x = -^29/2/10. 

o Case (142) A = i\: Every point of E belongs to S c except the point 
d = i. The canard cycle presented in the previous case is persistent for this 
case (remember that this canard cycle is born from the bifurcation of the 
S— graph of Case 122- So, it radius does not tend to zero when A tends to 
So the non generic tangential singularity d = i behaves itself like a weak 
attractor focus. See Figure [20l 

o Cases (152) h < A < M 3 and (I62) 
13i — 14i replacing L\ by %\ and L3 by M3 
which Ti collides with T2- 

o Cases (17 2 ) M 3 < A < M 2 , (18 2 ) A = M 2 , (19 2 ) 
A = (5 and (21 2 ) A > ft: Analogous to Cases 15i — 19i 
and L3 by M3. 

The bifurcation diagram is illustrated in Figure [2TJ 



A = M 3 : 
where M3 



Analogous to Cases 
is the limit value for 

M 2 < A < p, (20a) 
replacing L 2 by M 2 

□ 



8. Proof of Theorem 3 

Proof of Theorem 3. In Cases I3, 23 and 33 we assume that Y presents the 
behavior Y~ . In Cases 43, 53 and 63 we assume that Y presents the behavior 
Y°. In Cases 7 3 - 21 3 we assume that Y presents the behavior Y + . 

o Cases (1 3 ) d\ < e\, (2 3 ) d\ = e\, (3 3 ) d\ > e x , (4 3 ) d\ < si, (5 3 ) 
d\ = s\ and (63) d\ > s±: Analogous to Cases li, 2i, 3i, 4i, 5i and 61. 



ON THREE-PARAMETER FAMILIES OF FILIPPOV SYSTEMS - THE FOLD-SADDLE SINGULARITY 




FIGURE 21. Bifurcation Diagram of Theorems 2 and 3. 



In what follows we consider Mq, Mi, Mi and M3 as in the previous 
theorem. 

o Cases (7 3 ) A < -f3, (83) A = -p, (9 3 ) -/3 < A < M 0j (10 3 ) A = M 
and (II3) Mo < A < ii: Analogous to Cases 72 — II2 changing M\ by £1. 

o Case (123) A = zi. - Every point of 2/{d} belongs to E c . In a similar 
way as Case 132, we can construct sequences (ki)i & ^ and (Zj)ieN- Since d = i 
we have that ki — > d and 2, — >• d. So d is a non generic tangential singularity 
that behaves itself like an attractor. See Figure [221 



A = h h<X<Mi A = Mi 




FIGURE 22. Cases 12 3 , 13 3 and 14 3 . 



o Case (133) i\ < A < Mi: Analogous to Case 132 except that there is a 
change of stability on P = (P, 0), which is a S— repeller, and on T, which 
is an attractor canard cycle of kind I. This canard cycle is born from the 
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bifurcation of the non-generic tangential singularity of Case 12%. See Figure 

E2 

o Case (143) A = M±: Analogous to Case 122 except that occurs a change 
of stability on P = (P,0), which is a S— repeller. This fact generates a 
bifurcation like Hopf near P and there appears a hyperbolic attractor canard 
cycle Ti, of kind I, between P and the X— graph T2- See Figure [221 

o Cases (153) M\ < A < M% and (I63) A = M3: Analogous to Cases 
152 — I62, replacing i\ by M\. 

o Cases (17 3 ) M 3 < A < Af 2 , (18 3 ) A = M 2 , (19 3 ) M 2 < A < /3, (20 3 ) 
\ = (3 and (2I2) A > /3: Analogous to Cases 172 — 2I2. 

The bifurcation diagram is illustrated in Figure (2TJ □ 

9. Proof of Theorem 4 

Proof of Theorem 4- Since X has a unique S— fold point which is visible we 
conclude that canard cycles do not arise. 

In Cases I4, 24 and 34 we assume that Y presents the behavior Y~. In 
Cases 44, 54 and 64 we assume that Y presents the behavior Y . In these 
cases, when it is well defined, the direction function H assumes positive 
values. 

o Case (I4) d\ < e±: The points of S inside the interval (di,e±) belong 
to S c . The points on the left of d\ belong to T, s and the points on the right 
of ei belong to S e . See Figure l23l 

o Case (24) d\ = e±: Here S c = 0. The vector fields X and Y are 
linearly dependent on d± = e\ which is a tangential singularity. Moreover, 
it is an attractor for the trajectories of Z crossing T, s and a repeller for the 
trajectories of Z crossing £ e . See Figure l23l 

o Case (34) d\ > e±: The points of E inside the interval (e±,di) belong 
to S c . The points on the left of e\ belong to S s and the points on the right 
of d\ belong to E e . See Figure [23l 




o Case (44) d\ < s\: The points of S inside the interval (d\,si) belong 
to S c . The points on the left of d\ belong to T, s and the points on the right 
of si belong to S e . See Figure [24l 
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o Case (54) d\ = s±: Here £ c = and S is an attractor for the trajectories 
of Z crossing T, s and it is a repeller for the trajectories of Z crossing E e . 
See Figure [241 

o Case (64) d\ > s±: The points of E inside the interval (d\,si) belong 
to S c . The points on the left of si belong to T, s and the points on the right 
of d\ belong to S e . See Figure l24l 



A < 



4 4 




/ 


\ 



5.1 




A = 



A > 



64 




FIGURE 24. Cases 4 4 , 5 4 and 64. 



In Cases 74 — 134 we assume that Y presents the behavior Y + . 

o Cases (74) d± < hi, (84) d\ = /ii and (94) h\ < d\ < i\: The points of 
S inside the interval (di, ii) belong to S c . The points on the left of d\ belong 
to S s and the points on the right of i\ belong to S e . The direction function 
H assumes positive values on T, s and negative values in a neighborhood of 
ii. Moreover, H((3X/(-l + /3)) = and the S-repeller P = (/?A/(-l + /3), 0) 
is the unique pseudo equilibrium. See Figure [25j 

o Case (IO4) d\ = i\: Here S c = 0. The vector fields X and Y are 
linearly dependent on the tangential singularity d\ = i\. A straightforward 
calculation shows that H(z) = (1 — /3)/2 7^ for all z 6 E/{d}. So d± = i\ 
is an attractor for the trajectories of Z crossing T, s and a repeller for the 
trajectories of Z crossing S e . Moreover, A = {d} U dj U S e U {S} U S c U hd 
is a S— graph of kind III in such a way that each Q in its interior belongs 
to another S— graph of kind III passing through d. See Figure [ 



A< -P 




10 4 









FIGURE 25. Cases 7 4 - 10 4 . 
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o Cases (II4) i\ < d\ < j\, (124) d\ = j\ and (134) j\ < d\: The points of 
£ inside the interval (ii, di) belong to S c . The points on the left of ii belong 
to S s and the points on the right of d\ belong to S e . The direction function 
H assumes positive values on S e and negative values in a neighborhood of i±. 
Moreover, H((3X/(-l + $)) = and the S-attractor P = (/3A/(-l + /3),0) 
is the unique pseudo equilibrium. See Figure [26l 



0<A</3 A = /3 < A 




Figure 26. Cases 11 4 - 13 4 . 




FIGURE 27. Bifurcation Diagram of Theorems 4, 5 and 6. 

The bifurcation diagram is illustrated in Figure [271 d 

10. Proof of Theorem 5 

Proof of Theorem 5. The direction function if has a root Q = (OS 0) where 
(14) 

g= 2(a + l) (( ~ 1 + a)(1 " /3) " A(1 + a)+ 

+ «)(!-/?)- A(l + a)) 2 + 4/3(1 + a)(l + a + A(-l + a))). 
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Moreover, H assumes positive values on the right of Q and negative values 
on the left of Q. Note that when a —> —1 so Q — > — oo under the line 
{y = 0} and it occurs the configurations showed in Theorem 4. 

In Cases I5, 25 and 3s we assume that Y presents the behavior Y~ . In 
Cases 45, 55 and 65 we assume that Y presents the behavior Y°. In Cases 
7 5 — 135 we assume that Y presents the behavior Y + . 

o Cases (1 5 ) d\ < e 1; (2 5 ) d\ = e 1; (3 5 ) di > e\, (4 5 ) d\ < s 1; (5 5 ). 
d\ = si and (65) d\ > s±: Analogous to Cases I4, 24, 34, 44, 54 and 64 
respectively, except that here it appears the X— saddle Q on the left of d 
and e or S. See Figure 1281 



A < (1 + a)/3/(l - a) A = (1 + a)(3/(l - a) A > (1 + a)0/(l - a) 




2 5 




\/y\ 




FIGURE 28. Cases 1 5 , 2 5 and 3 5 



o Cases (7 5) d\ < h±, (85) di = h\, {%) h± < d± < i±: Analogous to 
Cases 74 — 94, except that here the X— saddle Q appears on the left of d\ 
and i\. So P = (p, 0) where 
(15) 



P 



■((-l + a)(l-/3)-A(l + a)+ 



2(a + 1) 

- V((-l + q)(1 - /3) - A(l + q)) 2 + 4/3(1 + «)(! + « + A(-l + a))). 



o Case (IO5) di = ii: Analogous to Case IO4, except that here appear 
the X— saddle Q on the left of d± = i±. 

o Cases (II5) h < d± < j±, (125) di = j\ and (13s) ji < d\: Analogous 
to Cases II4 — 134, except that here the X— saddle Q appears on the left of 
d\ and i\. 

The bifurcation diagram is illustrated in Figure [271 d 



11. Proof of Theorem 6 

Proof of Theorem 6. The direction function H has a root Q = (<?, 0) where 
q is given by (|14j) . Moreover, assumes positive values on the left of Q 
and negative values on the right of Q. Note that when a— 1 so Q — > 00 
under the line {y = 0} and the configurations shown in Theorem 4 occur. 
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In Cases 1q, 2q and 3q we assume that Y presents the behavior Y~ . In 
Cases Aq, 5q and 66 we assume that Y presents the behavior Y°. In Cases 
7q — 13q we assume that Y presents the behavior Y + . 

o Cases (1 6 ) d\ < e lt (2 6 ) di = e lt (3 6 ) d\ > ex, (4 6 ) d\ < s 1} (5 6 ) di = s\ 
and (6 6 ) d\ > s 1} (7 6 ) d\ < hi, (8 6 ) d\ = hi, {%) hi < di < i\, (10 6 ) 
di = h, (11 6 ) ii < di < ii, (12 6 ) di = ji anrf (13 6 ) ji < di: Analogous 
to Cases 1 5 , 2 5 , 3 5 , 4 5 , 5 5 , 65, 7 5 , 85, 9 5 , 10 5 , 11 5 , 12 5 and 13 5 respectively, 
except that here the S— saddle Q takes place on the right of d±, e±, s\ and 
i\ when these points appear. 

The bifurcation diagram is illustrated in Figure [271 d 
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